Journal of Radioanalytical Nuclear Chemistry manuscript No. 

(will be inserted by the editor) 



Investigation of production routes for the 6 Ho Auger-electron 
emitting radiolanthanide, a candidate for therapy 

F. Tarkanyi • F. Ditroi • A. Hermanne • S. Takacs • A.V. Ignatyuk 



o 

(N 
C 

o 



43 

Oh 
i 

-a 



c*3 
O 



> 
o 



Received: 2012-11-13 / Accepted: 2012-12-19 

Abstract The radiolanthanide 161 Ho (2.48 h) is a pro- 
mising Auger-electron emitter for internal radiother- 
apy that can be produced with particle accelerators. 
The excitation functions of the ™ at Dy(p,xn) 161 Ho and 
Tlat Dy(d,x) 161 Ho reactions were measured up to 40 and 
50 MeV respectively by using the stacked foil activa- 
tion method and -ray spectrometry. The experimental 
data were compared with results of the TALYS code 
available in the TENDL 2011 library I]. The main pa- 
rameters of different production routes are discussed. 

Keywords medical radioisotopes • therapeutical 
isotopes • proton and deuteron irradiation • 161 Ho • 

162m Ho 



1 Introduction 

The radiolanthanide 161 Ho is an Auger-electron emitter 
having also low energy photons in high abundance. It 
is very suitable for internal radiotherapy of small tu- 
mors because of the low energy electrons emitted. The 
short range of Auger electrons however requires that 
labeled compounds approach the cell nucleus. It is also 
interesting as low-energy narrow-band X-ray source for 
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internal irradiation [2,3,4,5,6 J. Different routes exist to 
produce 161 Ho with particle accelerators. One route is 
a- or 3 He- particle irradiation of 159 Tb relying on the 
159 Tb(a,2n) 161 Ho and 159 Tb( 3 He,n) 161 Ho reactions. An- 
other route is the irradiation of dysprosium targets us- 
ing protons via the 161 Dy(p,n), 162 Dy(p,2n) or deuterons 
via 160 Dy(d,n) and 161 Dy(d,2n) reactions. When using 
the 159 Tb(a,2n) 161 Ho, 162 Dy(p,2n) or 161 Dy(d,2n) re- 
actions processes, emission of one neutron also takes 
place resulting in simultaneous production of 162 Ho, 
which is a radionuclide impurity. This radionuclide has 
a short half-life ground state 1629 Ho (15 min, I" = 1 + ) 
and a longer-lived isomeric state (67.0 min, V = 6~). 
From the point of view of 161 Ho production, the con- 
tamination with the longer-lived excited state has some 
importance at the beginning after EOB, which will be 
reduced by the waiting time. The decay through inter- 
nal transition of 162m Ho is followed only by low energy 
low intensity 7-ray emission, but the 38% electron cap- 
ture decay results in strong high energy 7-lines (see Ta- 
ble 1). The cross-sections of the 159 Tb(a,2n) 161 Ho re- 
action was investigated by several authors (see in com- 
parison of production routes) . Although the basic cross- 
section data are still missing the proton induced reac- 
tion was also was used for practical production [5] . The 
deuteron induced reaction was not used yet for the pro- 
duction and no cross-section data were published. We 
decided to investigate the excitation functions of the 
proton and deuteron routes experimentally. Naturally 
occurring dysprosium is composed of 7 stable isotopes 
( 156 Dy - 0.06%, 158 Dy - 0.10%, 160 Dy - 2.34%, 161 Dy 



18.9%. 
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Dy - 25.5%, 163 Dy - 24.9% and 164 Dy 



28.2%). Taking into account that many of these sta- 
ble Dy isotopes have a larger mass number than 161 Ho, 
the direct experimental investigation of the 161 Dy(p,n), 
161 Dy(d,2n) or 160 Dy(d,n) reactions, and routes leading 
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Table 1 Main experimental parameters 



Reaction 


Mt Dy( P ,x) 


nat Dy(d,x) 








Incident particle 


Proton 


Deutcron 


Method 


Stacked foil 


Stacked foil 


Target and thickness 


naT Dy foil, 100.59 
fj.m 


nat Dy foil, 100.59 
fim 


Number of target foils 


15 


15 


Accelerator 


CGR 560 cyclotron 
of Vrijc Universiteit 
Brussels 


Cyclone 90 cyclotron 
of the Universit 
Catholique in City- 
placeLouvain la 
Neuvc (LLN) 


Primary energy 


36 McV 


50 McV 


Irradiation time 


71 min 


30 min 


Beam current 


61 nA 


120 nA 


Monitor reaction, 
[recommended values] 


nat Ti(p,x) 4 *V reac- 
tion Ell 


^'AlfdpO^Na reac- 
tion [T^] 


Monitor target and 
thickness 


nat Ti, 10.9 /Am 


nat Al, 26.96 fim 


detector 


HpGc 


HpGc 


-y-spectra measure- 
ments 


3 series 


3 series 


Cooling times 


1.5 h, 20 h, 80 h 


4h, 20h, 120h 




Energy (MeV) 



to possibly disturbing activation products, should re- 
quire highly enriched targets. A second possibility is to 
use natural targets and to use results of theoretical cal- 
culations to separate the contributions of the different 
target isotopes. This supposes an accurate predictivity 
of the calculations and a method to check the reliabil- 
ity has to be implemented. We adopted this approach 
by making measurement of production cross-section of 
161 Ho on Ilat Dy target with protons and deuterons and 
by comparing our experimental data with the predic- 
tions of the theoretical model codes. In case of good 
agreement we can then compare the different charged 
particle production routes using theoretical results val- 
idated by integral experiment. 



Fig. 1 The simultaneously measured monitor reactions for 
determination of proton beam energy and intensity 



Table 2 Main parameters of data evaluation (with 
references) 



7-spectrum evaluation 


Genie 2000, Forgamma 


I13I14I 


Determination of beam in- 
tensity 


Faraday cup (preliminary) 
Fitted monitor reaction 
(final) 


Il5l 


Decay data 


NUDAT 2.6 


1161 


Reaction Q-valucs 


Q-value calculator 


U71 


Determination of beam en- 


Andersen (preliminary) 


1181 


ergy 


Fitted monitor reaction 
(final) 


nr 


Uncertainty of energy 


cumulative effects of possi- 
ble uncertainties 


m 


Cross-sections 


Isotopic cross section 


[[20] 


Uncertainty of cross- 


Sum in quadrature of all 
individual contributions 


L±9] 


Yield 


Physical yield 


1201 



2 Experimental and data evaluation 

The general characteristics and procedures for irradia- 
tion, activity assessment and data evaluation (includ- 
ing estimation of uncertainties) were similar as in our 
many earlier works [TlfcTl^l lTCJlfTTj . The main experimen- 
tal parameters for the present study including the cho- 
sen monitor reactions |12j are summarized in Table 1. 
The main methods used in data evaluation and the used 
decay data p i3 p i4 | ll5 y i6 1 li7 1 [i8 l fT9 l l20] are collected in 
Table 2 and Table 3. The excitation function of simul- 
taneously measured proton and deuteron monitor re- 
actions and comparison with recommended values are 
shown in Fig. 1. 



3 Results 

3.1 Cross-sections 

The measured excitation functions for ™ at Dy(p,xn) 161 ' 162m - 
Ho and nat Dy(d,xn) 161 < 162m Ho are shown in Figs 2-3 
and 5-6 in comparison with the results of the model 



calculations. The numerical data important for further 
data evaluation are collected in Tables 4 and 5. In both 
cases the theoretical results reproduce exceptionally well 
the shape of the measured excitation functions, but an 
overestimation over the whole energy range is seen for 
the ™ at Dy(p,xn) 161 ' 162m Ho and nat Dy(d,xn) 161 Ho re- 
actions. In case of the ™ at Dy(d,xn) 162m Ho the agree- 
ment of the maximum is acceptable. The numerical 
values of theoretical results by a factor of 0.7 in case of 
a; Dy(p,xn) 161 ' 162m Ho and of 0.55 in case of x Dy(d,xn) 161 - 
Ho should be multiplied as a rough estimation. There 
is no normalization for :c Dy(d,xn) 162m Ho reaction. For 
further discussion we have normalized the theoretical 
cross-sections of the contributing reactions involved in 
161 Ho production with these factors. The comparison 
of the normalized TENDL 2011 cross-sections of the 
161 Dy(p,n) 161 Ho and the 162 Dy(p,2n) 161 Ho reactions 
and the 162 Dy(p,n) 162m Ho impurity reaction are shown 
in Fig. 4 and of the 160 Dy(d,n) 161 Ho, the 161 Dy(d,2n) 161 - 
Ho reactions and the 161 Dy(d,n) 162m Ho impurity reac- 
tion is presented in Fig 7. 
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Table 3 Decay characteristics of the 161 Ho and 162 Ho 
and Q-values of reactions for their productions) 



Nuclide 


Half- 
life 


B 7 (keV) 


M % ) 


Contributing 


Q-valuc 
(keV) 


^Ho 

e: 100 % 
7/2- 


2.48 h 


77.42 
103.05 

175.42 


1.9 

103.05 
0.43 


ltil Dy(p,„) 

162 Dy(p,2„) 

163 Dy(p,3„) 

164 Dy(p,4n) 

160 Dy(d,„) 

161 Dy(d,2n) 

162 Dy(d,3n) 

163 Dy(d,4n) 

164 Dy(d,5n) 


-1640.64 
-9837.63 
-16108.65 
-23766.77 

2589.18 

-3865.2 

-12062.2 

-18333.21 

-25991.33 


IT: 62 % 
e: 38 % 
6- 

105.87 kcV 


67.0 m 


184.99 
282.86 
937.17 
1220.04 


23.94 
10.27 
10.44 
23.7 


1 ' 2 D,(p l ») 

163 Dy(p,2„) 

164 Dy(p,3„) 

161 Dy(d,„) 
162 Dy(d,2n) 
163 Dy(d,3n) 
164 Dy(d,4„) 


-3027.91 
-9298.92 
-16957.05 

2944.51 
-5252.48 
-11523.49 
-19181.61 


iW 9Ho 

e: 100 % 
1 + 


15.0 m 


80.7 
1319.75 


8.0 
3.82 


lb2 Dy(p,n) 

163 Dy(p,2n) 

164 Dy(p,3„) 

161 Dy(d,n) 
162 Dy(d,2n) 
163 Dy(d,3n) 
164 Dy(d,4n) 


-2922.04 
-9193.05 
-16851.18 

3050.38 
-5146.61 
-11417.62 
-19075.74 




Proton energy (MeV} 

Fig. 2 Experimental cross-sections of the ™ at Dy(p,xn) 161 Ho 
reaction in comparison with the results of model calculations 
in TENDL 2011 




Proton energy (MeV} 

Fig. 3 Experimental cross-sections of the 

™ a *Dy(p,xn) 162rn Ho reaction in comparison with the 
results of model calculations in TENDL 2011 







161 Dy(p,n) 18!m Ho 
1s2 Dy(p,xn) 161 ' 162m Ho 




i 


162Dyft),n)162niHo 






162Dy(p i 2n) 161H0 




/ i 


ieiDy(p.n)1S1Ho 






\ 



Proton energy (MeV) 

Fig. 4 Experimental cross-sections of the ™ at Dy(d,xn) 161 Ho 
reaction in comparison with the results of model calculations 
in TENDL 2011 



Table 4 Measured cross-sections of the 
™ at Dy(p,x) 161162m Ho reactions 







l^m Ho 


E 


| ±Sa 


cr | ±8<T 


MeV 


mb 


mb 


35.3 


253.5 


34.8 


34.7 


5.7 


33.8 


246.2 


35.6 


46.9 


5.9 


32.3 


293.4 


39.0 


55.6 


7.3 


30.8 


285.5 


38.1 


82.6 


10.2 


29.2 


265.0 


35.2 


103.4 


11.5 


27.5 


214.8 


33.2 


113.6 


14.6 


25.7 


239.8 


32.2 


114.6 


13.1 


23.8 


239.2 


34.1 


118.6 


13.4 


21.9 


248.9 


28.1 


111.1 


12.3 


19.7 


226.1 


25.5 


93.0 


10.4 


17.4 


178.8 


20.6 


80.8 


9.0 


14.8 


156.1 


21.1 


73.0 


8.0 


11.9 


122.3 


14.1 


49.0 


5.4 


8.2 


39.2 


5.8 


15.3 


1.7 


4.3 


5.7 


1.7 


0.8 


0.1 



"»Dy(d,x) 1s1 Ho 




10 20 30 40 50 SO 70 80 90 100 

Deuteron energy (MeV) 



Fig. 5 Experimental cross-sections of the 

rlat Dy(d,xn) 162m Ho reaction in comparison with the 
results of model calculations in TENDL 2011 
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Table 5 Measured cross-sections 
nat Dy(d,x) 161162m Ho reactions 



of 



the 







162™ Ho 


E 


a | ±5cr 


er | ±8a 


MeV 


mb 


mb 


48. 6 


219.9 


26.6 


66.7 


9.0 


46.5 


222.2 


26.9 


75.2 


10.0 


44.4 


242.0 


28.1 


94.4 


11.2 


42.3 


295.2 


33.1 


93.9 


10.8 


40.0 


274. S 


31.3 


106.1 


12.5 


37.6 


301.0 


33.7 


143.2 


15.7 


35.1 


265. S 


29.8 


166.8 


18.4 


32.5 


277.7 


31.3 


192.5 


21.1 


29. S 


295.1 


33.1 


189.9 


20.8 


26. S 


250.9 


28.3 


172.8 


19.0 


23.6 


234.9 


26.3 


147.5 


16.3 


20.0 


181.0 


20.5 


126.1 


13.9 


16.0 


127.5 


14.7 


91.7 


10.1 


11.0 


45.0 


5.5 


36.6 


4.2 



/ 

/ 

/ 


\ ^Dyid.nJ^Ho 


/ 


\ 1B1 Dy(d,2n) 161 Ho 


/ 

/ 


\ 160 Dy(d,n) 161 Ho 
\ 

\ 


f 
/ 
/ 
/ 


\ 
\ 

\ 1SQDy(d,n)161 Ho 

\ 




\ 161 Dy(d.2n)161 Ho 

\ 




\ 161Dy{d.n)162mHo 

\ 




\ 




\ 




\ 




\ 




















Deuteron energy (MeV) 

Fig. 6 Comparison of the cross-sections of the 
161 Dy(p,n) 161 Ho and the 162 Dy(p,2n) 161 Ho reactions 
and the 162 Dy(p,n) 162 " l Ho impurity reaction in TENDL 
2011 



Deuteron energy (MeV) 

Fig. 7 Comparison of the cross-sections of the 
160 Dy(d,n) 161 Ho and the 161 Dy(d,2n) 161 Ho reactions 
and the 161 Dy(d,n) 162m Ho impurity reaction in TENDL 
2011 



Physical yields 










ratDvlMlSIHo 




— — natDY(p,x)152raHo 




natDv(d,x)161Ho 




natDy(d,xJ16ZmHo 




■ Stephens natDv(p,x]-161Ho 



3.2 Integral yields 



Particle energy (MeV) 

Fig. 8 Integral yields of the ™ at Dy(p,xn) 161 Ho, 
™ at Dy(p,xn) 162m Ho, ™ at Dy(d,xn) 162 " l Ho and 

™ Q *Dy(d,xn) 162m Ho reactions 



The integral yields calculated on the basis of the nor- 
malized TENDL 2011 cross-sections for the na *Dy(p,xn)- 
161 Ho, " at Dy(p,xn) 162m Ho, ™ Qt Dy(d,xn) 162m Ho and nat - 
Dy(d,xn) 162m Ho production reactions are shown in Fig. 
8. The calculated integral yield represents so called phys- 
ical yield i.e. yield obtained in a short irradiation [19 . 
The na *Dy(p,xn) 161 Ho yields are compared with the ex- 
perimental data of Stephens [5]. The value calculated 
by the Stephens' result for 11.6 MeV proton bombard- 
ment is significantly lower than our result, but it can 
be caused by the fact that the irradiation time was 
not published in that paper. Comparing the saturation 
activities, which are 1.8 GBq by Stephens (after 3.6 
hours) and 1.89 GBq in our measurement /calculation, 
the agreement can be considered as good. 



4 Comparison of production routes on different 
target materials 

The main parameters of the selected low and medium 
energy reactions that can lead to production of 161 Ho 
on different target materials are collected in Table 6. 
The excitation functions of the proton and deuteron 
routes are shown in Figs. 2-3 and Fig 4-5. Mukhcr- 
jee [21] and Bonesso [22] reported earlier experimental 
cross-section data on 159 Tb(a,2n) 161 Ho and Mukher- 
jee [H] and Singh [53] on total cross-section of and 
159 Tb(a,2n) 162 Ho. The 162m Ho/ 1629 Ho isomeric ratio 
was measured by Tulinov [21], Baskova [2S]- No exper- 
imental data were found for the 159 Tb( 3 He,n) 161 Ho re- 
action. The experimental data from literature and the 
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theoretical excitation functions of the 159 Tb(a,2n) 161 Ho 
and 159 Tb(a,n) 162m Ho reactions are shown in Fig. 9 
and 10 respectively. The theoretical excitation functions 
of the 159 Tb(a,2n) 161 Ho and 159 Tb(,n) 162m Ho and the 
159 Tb( 3 He,n) 161 Ho reactions are compared in Fig. 11. 
From the excitation functions of the above mentioned 
reactions the following conclusions can be drawn: 

— The production yields for the 162 Dy(p,2n) is the 
highest followed by the 161 Dy(d,2n), 159 Tb(a,2n) 
and 161 Dy(p,n) 

— No 162m Ho impurity is produced when using of 159 Tb- 
( 3 He,n), 161 Dy(p,n) and the 160 Dy(d,n) reactions. 
Among them the 161 Dy(p,n) reaction has the high- 
est cross-section ( max 260 mb) followed by the 
160 Dy(d,n) reaction( max 60 mb) and the less pro- 
ductive 159 Tb( 3 Hc,n) ( max 1 mb) 

— The element Tb is monoisotopic, relatively cheap 
and recovery is practically not necessary 

— In case of proton and deuteron induced reactions 
highly enriched 160, 161 or 162 Dy targets are re- 
quired 

— The production cross-sections of the 162m Ho from 
169 Tb(a,n), 161 Dy(d,n) and 162 Dy(d,n) are low 

— The impurity level depends on the selected energy 
range. The ratio of the saturation yields of the main 
reaction and of competing impurity reaction is shown 
in Fig. 12 as a function of energy. In the production 
energy range the ratio is lower than 3 

— The half-life of 162rn Ho is three times shorter, there- 
fore by using a short irradiation, the activity im- 
purity level will reach 3%. But by using irradiations 
lasting two half-life of 161 Ho and taking into account 
1 hour needed for the chemical separation and the 
labeling process the impurity level of 162m Ho will be 
reduced to 1% by decay. 



15s Tb(a,2n) 151 Ho 




10 15 20 25 30 35 40 45 50 55 60 

Particle energy (MeV) 

Fig. 9 Experimental and theoretical cross-sections of the 
159 Tb(a,2n) 161 Ho reaction 



40 




f~\ 159 Tb(a,n) 162 Ho 


30 










\ x 162HO total Mukherjee 1997 


20 




\ 162Ho total TENDL 201 1 






\ 162mHo TENDL 2011 






j \ \ a 1S2Ho total Singh 2001 


10 







10 15 20 25 30 35 40 45 50 55 60 

Alpha particle energy (MeV) 

Fig. 10 Experimental and theoretical cross-sections of the 
159 Tb(a,2n) 162m Ho reaction 



5 Summary and conclusions 

The principal aim of this work was an investigation of 
the production possibility of the radiotherapy related 
161 Ho. We present first experimental cross-sections for 
™ at Dy(p,xn) 161 Ho, ™ Qt Dy(p,xn) 162m Ho, rla *Dy(d,x) 161 Ho 
and " a *Dy(d,x) 162m Ho up to 40 and 50 MeV incident 
particle energies respectively. The TENDL 2011 the- 
oretical data predict well the shape of the excitation 
functions but overestimate the absolute values with a 
nearly constant factor in the whole energy range. The 
comparison of the different production routes shows 
that for production of 161 Ho with high radionuclide pu- 
rity the 161 Dy(p,n) 161 Ho, 162 Dy (p,2n) and 161 Dy(d,2n) 
reactions give the highest production yields. The 162m Ho 



radionuclide impurity level of the last two reactions 
however is significant. No enriched target material is 
necessary in case of 159 Tb(a,2n) (Tb is monoisotopic) 
but it requires accelerators having medium energy al- 
pha particles. The 159 Tb( 3 He,n) and 160 Dy(d,n) reac- 
tions have very low cross-sections, and accelerators dis- 
posing of 3 He beam are rare and the 3 He irradiation 
without recovery of 3 He gas is expensive. On the basis 
of the production yields, the impurity levels and the 
requirements of the medical application the 161 Dy(p,n) 
reaction is the production method of best choice. 
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of the HAS-FWO Vlaanderen (Hungary-Belgium) project. 
The authors acknowledge the support of the research project 
and of the respective institutions in providing the beam time 
and experimental facilities. 
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Table 6 Summary of the production parameters for selected reactions 



Reaction 


Q-value 


Impurity reac- 
tion 


Optimal 
energy range 
(MeV) 


lbl Ho 
thick tar- 
get yield 
(GBq/C) 


Impurity 
level 

(%) 


Optimal 
energy range 
at low impu- 
rity 
(MeV) 


lbi Ho 
thick tar- 
get yield 
(GBq/C) 


Impurity 
(%) 


I5y Tb(a,2n) lbi Ho 


-16053.93 


lba Tb(c<,n) Ibi! ™Ho 


35-19 


165 


6 


30-23 


107 


2.8 


lba Tb( 3 He,n) lbl Ho 


4523.7 




30-15 













Ibl Dy(p,n) Ibl Ho 


-1640.64 




15-8 


132 











Ib:! Dy(p,2n) lbl Ho 


-9837.63 


" 2 Dy( Pl n) ll;!m Ho 


30-12 


1459 


8.5 


22-15 


868 


4.6 


ibb Dy(d,n) ibl Ho 


2589.18 




15-7 


34 











Ibl Dy(d,2n) lbl Ho 


-3865.2 


lbI Dy(d,n) lbi!m Ho 


30-16 


1454 


1.8 


28-20 


955 


1.5 



15s Tb(a,2n) 1s1 Ho 

159Tb(3He,n)161Ho 
159Tb(a,2n)161Ho 
159Tb(a,n)162mHo 




Alpha- and ^He-particle energy (MeV) 

Fig. 11 Comparison of the cross-sections of the 
159 Tb(a,2n) 161 Ho and the 159 Tb( 3 He,n) 161 Ho reac- 
tions and the 159 Tb(a,n) 162m Ho impurity reaction (TENDL 
2011) 







1 

| "^Dy/^Dy yields 






1 162Cy(p,n)/162Dy(p,2n) 






161D/(d,n)/16IDy(d,2n) 






\ - 159Tb(a,n)/159Tl>;a,2n) 






\ 




V \ 


\ 




V V 


\ 















Particle energy (MeV) 



Fig. 12 12. The ratio of the saturation yield of the main 
reaction and the satellite impurity reaction 
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